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£\j , Abstract. Recently, Lembert, Gilson et al proposed a lucid and system- 



atic approach to obtain bilinear Backlund transformations and Lax pairs 
for constant-coefficient soliton equations based on the use of binary Bell 
polynomials. In this paper, we would like to further develop this method 
with new applications. We extend this method to systematically investigate 
complete integrability of nonisospectral and variable-coefficient equations. In 
addiction, a method is described for deriving infinite conservation laws of 
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nonlinear evolution equations based on the use of binary Bell polynomials. All 
conserved density and flux are given by explicit recursion formulas. By taking 
variable-coefficient KdV and KP equations as illustrative examples, their 
bilinear formulism, bilinear Backlund transformations, Lax pairs, Darboux 
covariant Lax pairs and conservation laws are obtained in a quick and natural 



manner. In conclusion, though the coefficient functions have influences on a 
variable-coefficient nonlinear equation, under certain constrains the equation 
turn out to be also completely integrable, which leads us to a canonical 
interpretation of their iV-soliton solutions in theory. 
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1. Introduction 

In many physical situations, it is often preferable to have an equation 
with variable-coefficients, which may allow us to describe real phenomena 
in physical and engineering fields. For example, variable-coefficient non- 
linear Schrodinger-typed ones, which describe such situations more re- 
alistically than their constant-coefficient counterparts, in plasma physics, 
arterial mechanics and long-distance optical communications PQ-|5]|. Many 
physical and mechanical situations governed by variable-coefficient KdV 
(vc-KdV) equation, e.g., the pulse wave propagation in blood vessels and 
dynamics in the circulatory system, matter waves and nonlinear atom 
optics enhanced by the observations of Bose-Einstein condensation in 
the weakly interacting atomic gases, the nonlinear excitations of a Bose 
gas of impenetrable bosons with longitudinal confinement, the nonlinear 
waves in types of rods [6]-|llj. In recent years, there has been consid- 
erable interest in the study of variable-coefficient nonlinear equations, 
such as vc-KdV, vc-KP, vc-Schrodinger, vc-Boussinesq and cylindrical 
KdV equations. Recent progress in the investigation of the complete in- 
tegrability and exact solutions for such equations via Painleve analysis, 
inverse scattering transformation, Hirota bilinear method and Darboux 
transformation has been reported, the details can be seen in reference 
and references therein [H]-[{21]. It is obvious that variable-coefficient 
equations are often more complicated and difficult to be solved than 
constant-coefficient ones. As well-known, investigation of integrability 
for a nonlinear equation can be regarded as a pre-test and the first step 
of its exact solvability. There are many significant properties, such as 
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Lax pairs, infinite conservation laws, infinite symmetries, Hamiltonian 
structure, Painleve test that can characterize integrability of nonlinear 
equations. This may pave the way for constructing their exact solu- 
tions explicitly in a future. But in contrast to constant-coefficient cases, 
very little of detail is known about complete integrability of variable- 
coefficient nonlinear equations. 

Among the direct algebraic methods applicable to nonlinear partial 
differential equations in soliton theory, there is on which has proved par- 
ticularly powerful: the bilinear method developed by Hirota [23, 125] . 
Once a nonlinear equation is written in bilinear forms by a dependent 
variable transformation, then multi-soliton solutions are usually obtained 
[29]-[j34]. The search for a Hirota representation of a given nonlinear 
equation is generally recognized as an important first step in the con- 
struction of multi-soliton solutions. Yet, the construction of such bilinear 
Backhand transformation is not as one would wish. It relies on a par- 
ticular skill in using appropriate exchange formulas which are connected 
with the linear representation of the system. Recently, Lembert, Gilson 
et al proposed an alternative procedure based on the use of Bell polyno- 
mials which enable one to obtain parameter families of bilinear Backhand 
transformation for soliton equations in a lucid and systematic way [36]- 
| I38| . The Bell polynomials are found to play an important role in the 
characterization of bilinearizable equations. As a consequence bilinear 
Backhand transformation with single field can be linearize into corre- 
sponding Lax pairs. Their method provides a shortest way to bilinear 
Backhand transformation and Lax pairs of nonlinear equations, which 
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establishes a deep relation between integrability of a nonlinear equation 
and the Bell polynomials. 

The problem that we consider in this paper is to further develop this 
method with new applications. We extend the binary Bell polynomi- 
als approach to a large class of nonisospectral and variable-coefficient 
equations, such as nonisospectral and variable-coefficient KdV and KP 
equations etc. One of the many remarkable properties that deemed 
to characterize soliton equations is existence of an infinite sequence of 
conservation laws. Here we propose a approach to construct infinite 
conservation laws of nonlinear evolution equations through decoupling 
binary Bell polynomials into a Riccati type equation and a divergence 
type equation. As illustrative examples, the bilinear representations, bi- 
linear Backlund transformations, Lax pais and infinite conservation laws 
of the vc-KdV and vc-KP equations are obtained in a quick and natural 
manner. The integrable constraint conditions on the variable-coefficient 
functions can be naturally found in the procedure of applying binary Bell 
polynomials. We can also find that though the coefficient functions have 
influences on a variable-coefficient equation, under certain constrains the 
equation still can admit many integrability properties which are similar 
to those of its standard constant-coefficient equation. The organization 
of this paper is as follows. In section 2, we briefly present necessary 
notations on multi-dimensional binary Bell polynomial that will be used 
in this paper. In the sections 3 and 4, we deal with integrability of non- 
isospectral vc-KdV equation and vc-KdV equation, respectively. We aim 
at integrability of nonisospectral vc-KP equation and vc-KP equation in 
the sections 5 and 6, respectively. 



2. Multi-dimensional binary Bell polynomials 
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The main tool used in this paper is a class of the Bell polynomials, 
named after E. T. Bell [35]. To make our presentation easy understand- 
ing and self-contained, we simply recall some necessary notations on the 
Bell polynomials, the details refer, for instance, to Lembert and Gilson's 
work PJ-EB]. 

Let / = /(#!,••■ ,x n ) be a C°° function with multi-variables, the 
following polynomials 

Y (f) = Y (f ) = e~ f d Ul ■ ■ ■ d ne e f 

1 mxi,— ,n e x£ \J J — 1 n u — ,nt\jrixi,— ,r e x e J u x e c 

is called multi-dimensional Bell polynomials, in which we denote that 
fnx 1: - , rt xt = ® r x\ ■ ■ ■ d x e f, ri = 0, • ■ ■ , ni; • • ■ ; r e = 0, • • • ,n e . 

For example, for the simplest case / = f(x), the associated one- 
dimensional Bell polynomials read 

Yi(f) = fx, Y 2 (f) = hx + fl Y 3 (f) = f 3x + 3fJ 2x + fl ■ ■ ■ . 

For / = f(x,t), the of associated two-dimensional Bell polynomials are 

Y x ,t{f) = fx.t + fxfti Y2 x ,t{f) = f2x,t + f2xft + %fx,tfx + fxfti 

Based on the use of above Bell polynomials, the multi-dimensional 
binary Bell polynomials can be defined as follows 

y ni x 1 ,-,n l x l (v,w) = Y nir . M (f)\ (v riXu ... m , n + Yrt is odd, 

fn x i>- > r t x t = \ 

[w nxir - ,r t x t , r l H ^ r £ is even > 

which inherit the easily recognizable partial structure of the Bell poly- 
nomials. 
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The first few lowest order binary Bell Polynomials are 

y x (v) = v x , y 2x {v, w) = w 2x + v%, y x ,t(v, w) = w xt + v x v t . 

(2.1) 

^3x = v 3x + 3v x w 2x + vl, ■ ■ ■ . 

The link between 3^-polynomials and the standard Hirota bilinear 
equation D™* • • • D^F • G can be given by an identity 

y niXlr . m (v = \nF/G,w = In FG) = (FG)~ l D™\ ■ ■ ■ D%F ■ G, (2.2) 

in which n\ + n 2 + • — h ri£ > 1. In the particular case when G = F, the 
formula (2.2) becomes 

F- 2 D n x l ■ ■ ■ D n x \F ■ F = y niXl ,.. M (0, q = 2 In F) 

( 0, m H h n £ is odd, (2.3) 

t^man,-^^), nH is even, 

in which the P-polynomials can be characterized by an equally recog- 
nizable even part partitional structure 

P 2x {q) = q 2x , Px,M) = Qxt, P^ x {q) = q^x + 3^, 
PqxW = q&x + 15q 2x q4 X + ^q 2x , 

The formulae (2.2) and (2.3) will prove particularly useful in connecting 
nonlinear equations with their corresponding bilinear equations. This 
means that once a nonlinear equation is expressible as a linear combina- 
tion of P-polynomials, then it can be transformed into a linear equation. 



7 

It follows that the binary Bell polynomials ^ niXl ..., w (f, w) can be 
separated into P-polynomials and y-polynomials 

(FG)~ 1 D™ 1 i ■ ■ ■ D n x l F ■ G = y ni xu-,ntxt(v, w)\ v= \ nF / GiW= \ nFG 

= y ni x u -,n e xA V y V + Q>)\v=]nF/G,q=2hiG 

= ' ' ' \\ ( r . ) Prixx,- ,r e xe(o)Y(ni-ri)xi,- ,(ne-re)x e ( v ) ■ 



niH Ym=even n=0 r^=0 i=l 



(2.5) 

The key property of the multi-dimensional Bell polynomials 

V Ul _ VVxi,--- ,n e xe (Q p>\ 

1 mxu- ,n e x t \u ) \v=inip — , , l z - u J 

implies that the binary Bell polynomials y niXl ,- ,nix e ( v > w) can still be 
linearized by means of the Hopf-Cole transformation v = \nip, that is, 
ip = F/G. The formulae (2.5) and (2.6) will then provide the shortest 
way to the associated Lax system of nonlinear equations. 

We start with construction of infinite conservation laws by virtue of 
binary Bell polynomials. We define a new auxiliary field variable 

V= (Q Xk -Qx k )/2, 

where q' and q are given hy q 1 = w + v, q = w — v and Xk is a appropriate 
variable chosen from xi, • • • ,X£. The two- filed condition 

C(q', q) = E(q) — E{q) = (2.7) 

can be regarded as the natural ansatz for a bilinear Backbend transfor- 
mation. By expressing the two-filed condition (2.7) in terms of binary 
Bell ^-polynomials and their derivatives, we expect that the resulting 
condition is then decoupled into a pair of constrains, i.e. often a Riccati 
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type equation with respect to Xk 



r) Xk + fin) = 0, 



(2.8) 



and a divergence-type equation 



(2.9) 



The recursion formulas of conversed density come from the equation 

(2.8) , the formulas of associated flux are obtained by using the equation 

(2.9) . It is often the case that the first few conservation laws of a non- 
linear equation have a physical interpretation. 

3. Nonisospectral variable coefficient KdV equation 

Consider nonisospectral vc-KdV equation |T3] 



where hi = hi(t), hi = hi(t) and h% = h^(t) are all arbitrary func- 
tions with respect to time variable t. The equation (3.1) includes some 
governing physical equations as special reduction, such as celebrated 
constant-coefficient KdV equation 

u t + 6uu x + u ix = 0, 

cylindrical KdV equation (hi = 1, hi = l/8t, /13 = 0) [T3] 




(3.1) 



u t + 6uu x + u 3x + — u x = 
and the vc-KdV equation (hi = 1, hi = Cq/4, /i3 



(3.2) 



7) 



u t + 6uu x + u 3l + ju + [(cq + -yx)u] x = 0, 



which describes the effect of relaxation inhomogeneous medium [15] . It 
can be observed that the equation (3.1) is invariant under Galiean trans- 
formation 

u — > u + X, x — > x + 6Xt, t — >• t. 

The inverse scattering transformation of the equation (3.1) was consid- 
ered by Chan and Li [H]. Lou and Ruan obtained infinite conservation 
laws M17M . Here we shall investigate the integrability of the equation 



(3.1) from bilinear representation, Backlund transformation, Lax pair, 
Darboux covariant Lax pair and infinite conservation laws. 

3.1. Bilinear representation 

In order to detect its existence of linearizable representation, we in- 
troduce a potential field q by setting 

u = c(t)q 2xj (3.3) 

with c = c(t) being free function to be the appropriate choice such that 
the equation (3.1) connect with P-polynomials. Substituting (3.3) into 
(3.1) and integrating with respect to x yields 

E(q) - q x t + hi(q± x + 3cql x ) +4h 2 q 2x - h 3 (q x + xq 2x ) +q x d t \nc = 0. (3.4) 

Comparing the second term of this equation together with the formula 
(2.4), we require c(t) = 1. The result equation is then cast into a com- 
bination form of P-polynomials 

E(q) = P xt (q) + ftiPtofa) + 4h 2 P 2x (q) - h 3 (xP 2x (q) + q x ) = 0. (3.5) 

Making a change of dependent variable 

q = 2 In F u = cq 2x = 2 (In F) 2x 
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and noting the property (2.3), the equation (3.5) gives the bilinear rep- 
resentation as follows 

(D x D t + hiDi + 4h 2 D 2 x - xh 3 D 2 x - hid x )F ■ F = 0, 

in which we have used the notation d x F ■ F = d x F 2 = 2FF X . This 
equation is easy to be solved for multi-soliton solutions by using Hirota's 
bilinear method. For example, the regular one-soliton like solution reads 

k 2 okx + lu 

u = — sech , 

2 2 ' 

where k = k(t) and lu = u(t) are two functions about t, given by 

k(t) = aef h3dt , u(t) = - J (hk 3 + 4h 2 k)dt. 

The multi-soliton solution are omitted here since exactly solving the 
equation (3.1) is not our main purpose in this paper. 

3.2. Backlund transformation and Lax pair 

Next, we search for the bilinear Backlund transformation and Lax 
pair of the vc-KdV equation (3.1). Let q and q' be two different solutions 
of the equation (3.4), respectively, we associate the two-field condition 

E(q') - E{q) = (q' - q) xt + h^q 1 - q) 4x + 3(q' + q) 2x (q' - q) 2x ] 

(3.5) 

+ 4h 2 (q f - q) 2x - h 3 [(q' - q) x + x(q' - q) 2x ] = 0. 

This two-field condition can be regarded as the natural ansatz for a 
bilinear Backlund transformation and may produce the required trans- 
formation under appropriate additional constraints. 

To find such constraints, we introduce two new variables 

v = (q'-q)/2, w = (q' + q)/2, (3.6) 
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and rewrite the condition (3.5) into the form 
E(q') - E(q) = v xt + hi{v Ax + 6v 2x w 2x ) + 4h 2 v 2x - h 3 (v x + xv 2x ) 

(3.7) 

= d x [y t (v) + hM x (v, w)} + R(v, w) = 0, 
with 

R(v,w) = 3hiWionskmn[y 2x (v,w),y x {v)} + d x [4h 2 y x {v) - xh 3 y x (v)}. 

In order to decouple the two-field condition (3.7) into a pair of con- 
straints, we impose such a constraint which enable us to express R(v, w) 
as the x-derivative of a combination of 3^-polynomials. The simplest 
possible choice of such constraint may be 

y 2x (v,w) + ay x (v) = \, (3.8) 

where a and A are arbitrary parameters. On account of the equation 
(3.8), then R(v,w) can be rewritten in the form 

r(v, w) = d^hxy^v) + 4h 2 y x (v) - xh 3 y x (v)]. (3.9) 

Then from (3.7)-(3.9), we deduce a coupled system of ^-polynomials 

y2x(v, w) + ay x (v) - A = 0, 

(3.10) 

d x y t (y) + d^hM^v, w) + (3/nA + 4h 2 - xh 3 )y x (y)] = 0. 

where prefer the second equation in the conserved form without inte- 
gration with respect to x, which is useful to construct conservation laws 
later. By application of the identity (2.2), the system (3.10) immediately 
leads to the bilinear Backlund transformation 
(D 2 X + aD x - X)F • G = 0, 

[D t + hiDl + (3/iiA + Ah 2 - xh 3 )D x + (3}F ■ G = 0, 
where (3 is a arbitrary parameter. 
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By transformation v = Ini/j, it follows from the formulae (2.5) and 
(2.6) that 

y x (v) = lf>x/i> } y2x(v, W) = q 2x + ^2zM 

y3x(v, w) = 2>q 2x ^ x /^ + ipzx/ip, y t {v) = ip t /i>, 
on account of which, the system (3.10) is then linearized into a Lax pair 
with double parameters about A and (3 

L 1 i/j=(d 2 x + ad x + q 2x )i/j = Xij, X t = 2h 3 X, (3.11) 

(d t + L 2 )^ [d t + h x dl + 3/^2* + X)d x + (4/i 2 - xhz)d x ]i/> (3.12) 

or equivalently replacing q 2x by u, 

^ 2a; + ai/; x + (u- X)i/j = 0, A< = 2h 3 X, 

ipt + [/ii(2u + 4A + a 2 ) + 4/& 2 - 2^3]^ + W ~ h\u x - ah x {u - X)]ip = 0. 
Starting from this Lax pair, the Darboux transformation and soliton-like 
solutions of the vc-KdV equation (3.1) can be established [14]. It is easy 
to check that the integrability condition 

[Li - A, <9 t + L 2 ]V = 

is satisfied if u is a solution of the vc-KdV equation (3.1) and nonisospec- 
tral condition X t = 2h%X holds. 

3.4. Infinite conservation laws 

Finally, we show how to derive the infinite conservation laws for vc- 
KdV equation (3.1) based on the use of the binary Bell polynomials. The 
conservation laws actually have been hinted in the two-filed constraint 
system (3.10), which can be rewritten in the conserved form 
y2 X (v,w) + ay x (v)-X = 0, 



(3.13) 

d t y x (v) + d^hM^v, w) + (3M + 4/i 2 - xh 3 )y x (v)} = 0. 



by applying the relation d x y t (v) = d t y x (v) = v 
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By introducing a new potential function 

V = (q' x - fe)/2, 
it follows from the relation (3.6) that 

v x = r}, w x = q x + ij. (3.14) 
Substituting (3.14) into (3.13), we get a Riccati-type equation 

T] x + r] 2 + q 2x = X = e 2 , (3.15) 
and a divergence-type equation 

m + d^hrkx + 6/11(77 + e)e 2 - 2^(77 + ef + (4/i 2 - 2^3) (?7 + e)] = 0, 

(3.16) 

where we have used the equation (3.15) to get the equation (3.16) and 
set A = e 2 . 

To proceed, inserting the expansion 

00 

r J = e + Y, I ^^---) £ ~ n ^ ( 3 - 17 ) 

n=l 

into the equation (3.15) and equating the coefficients for power of e, we 
then obtain the recursion relations for I n 

T 1 T 1 1 

h = ~Px = ~-U, h = ~P2x = ~Ux, 

1 » _ (3-18) 

In+l = ~~^Jn,x + / J Ikln-k), TL = 2,3, • ■ ■ , 

k=l 

By applying the nonisospectral condition 

X t = 2h 3 X => e t = h 3 e, 
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then substituting (3.17) into (3.16) yields 



oo 



7 ^~ n + d - 



n=l 



oo oo oo 



n=l n=l n=l 

= 0, 



oo oo 



+ (4/l 2 - 2^3) XI /n£ " ~ ^ X n ^ ^ 

n=l n=l 

which provides us infinite consequence of conservation laws 

In,t + F n , x = 0, n= 1,2,-.. . (3.19) 

In the equation (3.19), the conversed densities /4 s are given by formula 
(3.15) and the fluxes F' n s are given by recursion formulas explicitly 

Fi = -\ [hi(u 2x + 3w 2 ) + Ah 2 u - h 3 (xu + d^u)] , 

1 

^2 = I [^l(«3a: + §UU X ) + 4/^^ ~ ^(^ + 2u)\ , 

n 

F n = - 6/11 X hh+i-k - 2/ii X /j/y-ifc + (4/i 2 - ^3)4 

fc=l i+j+k=n 

+ nhsd^In, n = 3, 4, • • • . 

(3.20) 

We present recursion formulas for generating an infinite sequence of con- 
servation laws for each equation, the first few conserved density and as- 
sociated flux are explicit. The first equation of conservation law equation 
(3.19) is exactly the vc-KdV equation (3.1). The expressions (3.20) indi- 
cate that the fluxes F^s of the vc-KdV equation are not local, which are 
different from those of standard const ant- coefficient KdV equation. In 
conclusion, the vc-KdV equation (3.1) is complete integrable in the sense 
that it admits bilinear Backhand transformation, Lax pair and infinite 
conservation laws. 



4. Generalized variable-coefficient KdV equation 
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A more general example, we consider vc-KdV equation 

u t + hiu 3x + h 2 uu x + h s u x + h A u = 0, (4.1) 

where hj = hj(t), j = 1,2,3,4 are all arbitrary functions with respect 
to time variable t. Special cases of the equation (4.1) include cylindrical 
equation [22J-[21] 

u t + f(t)uu x + g(t)u 3x = 



and other special variable-coefficient equation J23, [2S] 



u t + at n uu x + bt m u ix = 0. 

Recently, Zhang et al obtained the bilinear form, Backhand transforma- 
tion and exact solutions for the equation (4.1) under the constrain [20] 

h x =coh 2 e-f h4dt . (4.2) 

Here we construct bilinear representation, Backlund transformation, Lax 
pair and conservation laws of the equation (4.1) based on the sue of 
binary Bell polynomials technique, which will be seen to be a natural way 
to find such a constraint (4.2). We find that the bilinear representation 
of the equation (4.1) existence without need of any constraint. The 
constraint only need it for construction of the Backlund transformation, 
Lax pair and conservation laws. 

4.1. Bilinear representation 

As before, we introduce a field q by setting 

u = c(t)q 2x , (4.3) 
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in which c = c(t) is free function to be determined. Substituting (4.3) 
into (4.1) and integrating with respect to x yields 

E(q) - q xt + htfte + ^h 2 cq\ x + h 3 q 2x + (h 4 + d t In c)q x = 0. (4.4) 

which can be cast into a combination form of P-polynomials by using 
the formula (2.4) 



E{q) = P xt {q) + hP^q) + h 3 P 2x (q) + (h 4 + d t In hh^q* = 0, (4.5) 



and using the property (2.3), then the equation (4.5) implies the bilinear 
form for the vc-KdV equation (4.1) as follows 



[D x D t + hiDi + h 3 D 2 x + (h A + dtlnhh^d^F ■ F = 0, (4.6) 



which is obviously more general than that obtained in [2D] , since we have 
no any constraint on the hi : h 2: and h±. Starting the bilinear equation 
(4.6), we can get multi-soliton solutions to the vc-KdV equation (4.1). 
For example, one-soliton solution takes the form 



if one chooses the function c(t) = 6h\h 2 ■ 
By transformation 



q = 2lnF <^> u = c(t)q 2x = Uh^ 1 (In F) 2x 



1 cy 

u = 6h\h 2 k sech 



2 kx + uj(t) 



2 



where k is a constant and cu(t) given by 





4.2. Backlund transformation and Lax pair 

In order to obtain the bilinear Backlund transformation and Lax pairs 
of the equation (4.1), let q, q' be two solutions of the equation (4.4) and 
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consider the associated two-field condition 

E{q) - E{q) = (q' - q) xt + h^q' - q) Ax + 3(q' + q) 2x (q' - q) 2x ) 

(4.7) 

+ h 3 (q' -q) 2x + {h A + d t In h x h 2 l ){q' - q) x ] = 0, 
which may produce the required bilinear Backhand transformation under 
an appropriate additional constraint. By introducing variables 

v = (q > -q)/2, w = (q > + q)/2 (4.8) 

we can rewrite the condition (4.7) as the form 
E{q') - E{q) = v xt + hi(v 4x + 6v 2x w 2x ) + h 3 v 2x + (h 4 + d t In hih 2 l )v x 

= d x [y t (v) + hiyte(v, w)} + R(v, w) = 0, 

(4.9) 

with 

R(v,w) = 3/iiWronskian[3^2x(^, w),y x (v)] + h 3 v 2x + (h 4 + d t \nhih 2 x )v x . 

In order to express R(v,w) as the x-derivative of a linear combination 
of 3^-polynomials, we choose a constraint 

y 2x {v,w) + ay x (v) = \ (4.10) 

where a and A are arbitrary parameters. Direct calculation gives 

R(v, w) = 3hi\v 2x + h 3 v 2x + (h 4 + d t In hih 2 l )v x , 

which can be written as ^-derivative of ^-polynomials 

R(v,w) = d^Xy^v) + h 3 y x (v) + (h A + ^lnMa >]. (4.11) 

From (4.10)-(4.11), we infer that 
y2x(v, w) + ay x (v) - A = 0, 

d^tM + d x [h x y 3x (v, w) + (3/nA + h 3 )y x (v) + (/i 4 + $ In h x h 2 l )v] = 0, 

(4.12) 
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which can be cast into a bilinear Backhand transformation by using the 
property (2.2) 

(D 2 x + aD x -\)F-G = 0, 

(4.13) 

[D t + h x Dl + (3/nA + h 3 )D x + p\F ■ G = 0, 
if we set the constraint 

h A + ^ln/ii/i^ 1 = =^ hi = c e-f h4dt h 2 , 

with cq being a arbitrary parameter. Without loss of generality, taking 
Co = 1/6, then 

c(t) = 6h 1 /h 2 = e-f hidt . 

As a = (3 = 0, the Backhand transformation (4.13) reduces to the one 
obtained in [20]. 

Making use of the Hopf-Cole transformation v = Inip and the formula 
(2.5)-(2.6), then the system (4.10) can be linearized into a Lax pair 

Liip = {d 2 x + ad x + q 2x )^ = (4.14) 

(d t + L 2 )tP = [d t + h x dl + (3hiq 2x + 3Xhi + h 3 )d x ]ijj = 0, (4.15) 
or equivalently, 

^2, + a^ + M^-A)^ = 0, 

^ + [/ii(2we^ Mi + 4A + a 2 ) + h 3 ]ip x + [(3 + aX - h x {u x ^ Kdt - aue$ Kdt )]iP = 0. 
This Lax pair can be used to construct Darboux transformation, inverse 
scattering transformation for soliton solutions. It is easy to check that 
the integrability condition 

[L x - A,d, + L 2 ]^ = 

is satisfied if u is a solution of the vc-KdV equation (4.1). 
4.3. Darboux covariant Lax pair 
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Let us go back to the vc-KdV equation (4.1) and the associated Lax 
pair (4. 14)- (4. 15). Assume that is a solution eigenvalue equation (4.14) 
(taking a = for simplicity). It is well-known that the gauge transfor- 
mation 

T = ( j)d x (j)- 1 = d x -(j, a = d x ]n(f> (4.16) 
map the operator L\ = d% + q 2x onto a similar operator: 

T^q) - \)T- 1 = L^q) - \ 

which satisfies the covariance condition 

Li{q) = Li(q = q + Aq), with Ag = 21n0. 

But it can verified that similar property does not hold for the evolu- 
tion equation (4.15). Next step is to find another third order operator 
L2,cov(q) with appropriate coefficients, such that dt J rL 2jC0Y (q) be mapped, 
by gauge transformation (4.16), onto a similar operator d t + L 2 ^ cov (q) 
which satisfies the covariance condition 

L 2 , CO y(q) = L 2 ,cov(q = q + Aq). 

Suppose that is a solution of the following Lax pair 

L\(J) = \<f>, 

( 4 -!7) 

(d t + ^2,cov)0 = 0, L 2 , CO v = + hd x + b 2 , 
where b\ and b 2 are functions to be determined. It suffice that we require 
that the transformation T map the operator d t + L 2fiOV onto the similar 
one 

T(d t + La.oovjT" 1 = dt + L 2 , cov , L 2 , cov = 4M* + hd x + b 2 , (4.18) 
where b\ and b 2 satisfy the covariant condition 

bj = bj(q) + Abj = bj(q + Aq), j = 1, 2. (4.19) 
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It follows from (4.16) and (4,18) that 

A&i = b l -bi = 12hKT x , Ab 2 = b 2 -b 2 = h jX + a Ah + \2h x o 2x , (4.20) 
and a satisfies 

a t + 4hia lx + b 2jX + Ab 2 a + b x o x = 0. (4.21) 

According to (4.19), it remains to determine b\ and b 2 in the form of 
polynomial expressions in terms of derivatives of q 

bj = Fj(q, q X: q 2x , q 3x , •••), j = 1, 2 

such that 

AFj = Fj(q + Aq, q x + Aq x , q 2x + Aq 2x , ■■■) - Fj(q, q x , q 2x , •••) = Abj, 

(4.22) 

with Aq rx = 2(hLq) rx , r = 1, 2, • • • , the Abj being determined by the 
relations (4.19). 

Expanding the left hand of the equation (4.22) , we obtain 

A6i = AFi = F u Aq + F hqx Aq x + F hq2x Aq 2x + ■■■ = \2h x o x = <Sh x Aq 2x , 

which implies that we can choose 

h = Fi(fe) = Qhiq 2x + ci(t), (4.23) 

with C\(t) being arbitrary function about t. 

From the eigenvalue equation in (4.17), we can find the following 
relation 

Q3x = ~o 2x - 2oo x . (4.24) 
Substituting (4.23) and (4.24) into (4.19) leads to 

Ab 2 = I2hiqs x + V2h\aa x + \2h x o 2x = Qhia 2x = S^Aq^. 
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The second condition 

AF 2 = F 2 , q Aq + F 2 , qx Aq x + F 2 , q2x Aq 2x + F 2mx Aq ix + ' ' ' = A6 2 , 

can be satisfied, if one chooses 

&2 = F 2 (q, q x q 2x , q 3x ) = ?>h x q ix + c 2 (t), 

in which c 2 (t) is arbitrary constant. 

Setting c\{t) = /i3, C2(t) = 0, we find the following Darboux covariant 
evolution equation 

(d t + ^2,cov)0 = 0, L 2;CO v = 4hidl + + h 3 )d x + 

which is in agreement with the equation (4.21). Moreover, the relation 
between the operator L 2jCOV and the operator L 2 is given by 

L 2jC ov = L 2 + 3h\d x (Li — A). 

The integrability condition of the Darboux covariant Lax pair (4.15) 
precisely give rise to the equation (4.1) in Lax representation 

[d t + ^2,cov, Li] = -[q xt + hiq Ax + 2>hiq\ x + h 3 q 2x + h A q x ] x . 

The higher operators can be obtained in a similar way step by step 

L k , C oM = + hd k x - 2 + • • • + b p , k = 3, 4, • • • 

which are Darboux covariant with respect to Li, so as to produce higher 
order members of the vc-KdV hierarchy. 

4.4. Infinite conservation laws 

Finally, we construct the conservation laws of vc-KdV equation. The 
second equation of (4.12) has been conserved form due to the relation 
d x yt(v) = d t y x (v) = v xt . By introducing a new potential function 
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it follows from the relation (4.8) that 

v x = 77, w x = q x + 77. (4.25) 

Substituting (4.25) into (4.12), we get a Riccati type equation 

Vx + rf + q 2x = A = e\ (4.26) 

and a divergence type equation 

Vt + dxih^x + 6/11(77 + e)e 2 - 2/11(77 + e) 3 + /i 3 (77 + e)] = 0, (4.27) 

where we have used the equation (4.26) to get the equation (4.27) and 
set A = e 2 . 

Substituting the expansion 

00 

r) = e + Y,U<l,<h»- '•)*-"■ (4-28) 

n=l 

into the equation (4.26) and equating the coefficients of e -1 , the con- 
served densities are explicitly obtained by recursion relations 

h = -Px = -\e Sludt u, h = -/ hidt u x , 

1 » _ (4-29) 

I-n+l = ~ ^ (Jn,x + / J Ikln-k), n = 2, 3, • • • , 

fc=l 

In addition, substituting (4.28) into (4.27) leads to 



00 



n=l 



n=l n=l n=l 

which provides us infinite conservation laws 

I n ,t + F n , x = 0, n= 1,2,-.. . (4.30) 
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In the equation (4.30), the conversed densities I' n s are given by recursion 
formulas (4.29), and the fluxes F' n s are given by 

Fi = -^ h4dt (h lU2x + 3h 2 u 2 + h 3 u), 

F 2 = ^ef h * dt (h lU3x + 6h 2 ef h * dt uu x + h 3 u x ), 

n 

F n = hj ni2x - 6hi ^2 hln+i-k - 2/&i ^2 IJjh + h 3 I n , n = 3, 4, • • • . 

fc=l i+j+k=n 

The first equation of (4.30) is exactly the vc-KdV equation (4.1). We 
see that above fluxes F' n s can be obtained from solution u by algebraic 
and differential manipulation, thus they are local. Taking the bound- 
ary condition of u into account, the conservation equation (4.30) implies 
that {/ n , n = 1,2, ••• ,} constitute infinite conserved densities of the 
vc-KdV equation (4.1). In conclusion, the vc-KdV equation (4.1) is 
complete integrable under the constraint h 2 = §h\e$ h7dt in the sense 
that it admits bilinear Backhand transformation, Lax pair and infinite 
conservation laws. 

5. Nonisospectral variable-coefficient KP equation 



Consider nonisospectral and vc-KP equation [16 



u t + h\{u 3x + 6uu x + 3a 2 d x l u yy ) + h 2 (u x — axu y — 2ad x 1 u y ) 

(5-1) 

— h 3 (xu x + 2u + 2yu y ) = 0, 

where hi = hi(i), h 2 = h 2 (t) and h 3 = h 3 (t) are all arbitrary functions 
with respect to time variable t. The equation (5.1) reduces to the vc- 
KdV equation (3.1) when u = u(x,t) is independent of the variable y. 
In the case h\ = 1, h 2 = h 3 = 0, it reduce to standard KP equation. 
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5.1. Bilinear representation 

By introducing a potential field q 

u = c{t)q 2x , (5.2) 

with c = c(t) is a free function about t to be determined, the resulting 
equation (5.1) for q (integrating with respect to x) reads 

E(q) = q x t + hi(q Ax + Zcq\ x + Sa 2 q 2y ) + h 2 (q 2x - xaq xy - aq y ) 

(5.3) 

- h 3 {q x + xq 2x + 2yq xy ) + q x d t In c = 0, 
which can be expressed in the form of P-polynomials 
E(q) = P xt {q) + hi[P± x (q) + 3a 2 P 2y (q)] + h 2 [P 2x (q) - axP xy (q) - aq y ] 

- h 3 [xP 2x (q) + 2yP xy (q) + q x ] = 0, 

(5.4) 

if one chooses c(t) = 1 and use the formula (2.4). 
By application of the variable transformation 

g = 21nF u = 2(\nF) 2x 

and using the property (2.3), then the equation (5.4) give the bilinear 
form for the vc-KP equation (5.1) as follows 

[D x D t +h 1 (Di+3a 2 Dl)+h 2 (Dl-axD x D y -ad y )-h 3 (xDl+2yD x D y +d x ^ = 0, 

starting from which, we can get multi-soliton solutions to the equation 

(5.1). For example, one-soliton solution takes the form 

(k + s) 2 2 ^ + C + lna; 
u = sech , 

in which £ = kx — k 2 y/a, £ = sx + s 2 y/a, while k = k(t), s = s(t) and 
lu = Lu(t) are all functions about t, satisfying 
k t = h 3 k — bk , St = /13s + 6s , 



u(t) = exp( y [h 3 (k 3 + s 3 ) - /i 2 (/c + s)]dt). 
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5.2. Backlund transformation and Lax pair 

In following we consider bilinear Backlund transformation and Lax 
pairs of the equation (5.1). Assume q and q' are two solutions of the 
equation (5.3), we consider the corresponding two-field condition 

e (q') ~ e (q) = id ~ Q)xt + hi[(q ~ q)Ax + 3(g' + q) 2x (q ~ q)2x + 3a 2 (q' - 
+ h 2 [(q - q) 2x ~ ax(q' ~ q) xy ~ a{<i ~ q) y ] ~ h 3 [x(q - q) 2x 
+ 2y(q' - q) xy + {q f - q) x ] = 0, 

(5.5) 

which may produce the required bilinear Backlund transformation under 
an appropriate additional constraint. 
By change of the variables 

v = (q'-q)/2, w = (q' + q)/2 (5.6) 

we rewrite the condition (5.5) in the form 

E(q') - E(q) = v xt + hi(v4x + 6v 2 xW2x + 3a 2 v 2y ) + h 2 (v 2x - axv xy - av y ) 

- h 3 {xv 2x + 2yv xy + v x ) = d x [y t (v) + h r y ix {v, w)} + R(v, w) = 0, 

(5.7) 

with 

R(v, w) = 3/iiWronskian[3^2x(^, w),y x (v)] + 3hia 2 v 2y + h 2 (v 2x — axv xy — 

- hz(xv 2x + 2yv xy + v x ). 

In order to express R(v,w) as the x-derivative of a linear combination 
of 3^-polynomials, we choose the constraint 



y2x(v,w) + ay y (v) = A, 



(5.8) 
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on account of which 
R(v,w) = d x {hi[3\y x (v) - 3ay xy {v,w)} + h 2 [y x {y) - axy y {v)} 

(5.9) 

-h 3 [xy x (v) + 2yy y (v)}}. 

Combining relations (5.7)-(5.9), we then obtain a pair of constraints 
in Bell polynomial form 

y2x(v,w) + ay y (v)-\ = 0, 

d t y x (v) + d^h^y^iv, w) - 3ay xy {v, w) + Z\y x (v)] + h 2 [y x {v) - axy y {v)} 
-h 3 [xy x (v)+2yy y (v)]} = 0, 

(5.10) 

which result in bilinear Backlund transformation 
(D 2 x + aD y - X)F-G = 0, 

[D t + h x (Dl - 3aD x D y + 3XD X ) + h 2 (D x - axD y ) (5.11) 

- h 3 (xD x + 2yD y ) + 0]F • G = 

by applying the property (2.2). 

It only remains to linearize the expression (5.10) with the formulae 
(2.5) and (2.6). By using the Hopf-Cole transformation v = ln^, the 
bilinear Backlund transformation (5.10) can be linearized into a Lax pair 

(d y + Li)ip = ip y + a~ l ^ 2x + a' 1 {q2x ~ A)V> = 0, X t = 2h 3 \ 

(d t + L 2 )ip = ip t + ^ipte + [xh 2 + 2ya~ 1 h 3 )ip 2x + {6hiq 2x + h 2 - xh 3 )ip x 

+ [3hiq 3x - 3ahiq xy + (q 2x - \){xh 2 + 2ya~ 1 h 3 )]ip = 0, 

which can be used to construct Darboux transformation, inverse scatter- 
ing transformation for getting soliton solutions. It is easy to check that 
the integrability condition 

[d y + L h d t + L 2 ]ij = 
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is satisfied if u is a solution of the vc-KP equation (5.1) and nonisospec- 
tral condition X t = 2h 3 X holds. 

5.3. Infinite conservation laws 

Finally, we precede to construct the conservation laws of vc-KP equa- 
tion. For this purpose, we decompose the two-filed condition (5.7) into 
x- and ^/-derivative of a linear combination of ^-polynomials. Let us go 
back to the two-field condition (5.7) and consider the following decom- 
position 

R{v, w) = d x [3hi\v x - 3ahiv x v y + h 2 (v x - axv y ) - h 3 (xv x + 2yv y )] 
+ dy(3hia 2 v y + 3h\avl). 

(5.11) 

It follows from (5.7), (5.8) and (5.11) that 
y2x{v, w) + ay y (v) - A = 0, 

iy x {v)) t + d x {h x [y 3x (v, w) - 3ay y {v, w) + 3Xy x (v)] + h 2 [y x (v) - axy y {v)} 
- h 3 [xy x (v) + 2yy y (v)]} + d y [3hia 2 y y (v) + 3h ia y%(v)] = 0, 

(5.12) 

which is slightly different from (5.10) and can produce desired conserva- 
tion laws. By introducing a new potential function 

v = {q x - &)/2, 

and it follows from the relation (5.6) that 

v x = , n, w x = q x + rj. (5.13) 

Substituting (5.13) into (5.12), we decompose the two-field condition 
(5.7) into a Riccati type equation 

r] x + rj 2 + ad-\ + q 2x -e 2 = 0, (5.14) 
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and a divergence type equation 



T} t + d x [hi(r]2x + 6??£ -2t/ - 6arjd x rj y ) + h 2 (ri - xad x r) y ) 

(5.15) 

- /i 3 (x7/ + 2?/d x 1 t/j / )] + 5 y (3/iia 2 5 x ^ + Shiarj 2 ) = 0, 



where we have used the equation (5.14) to get the equation (5.15) and 
set A = e 2 . 

Substituting the expansion 



oo 

r/ = £ + ^/nM*,---K n - (5-16) 

n=l 

Inserting the equation (5.16) into the equation (5.14), equation the co- 
efficients for power of e, then we have the recursion relations for I n 



h = -\<hx = ~\ u i h = ^(u2x + ad x 1 u y ), h = ~(u3x + u 2 + au y ), 



1 

4+i = -~(In,x + ^2 IkIn ~ k + udx 1 ^), n = 3, 4, 



2 

fc=i 



(5.17) 

Again substituting (5.16) into (5.15) and noting the nonisospectral con- 
dition 



X t = 2h 3 X e t = h 3 e, 



we then obtain the following infinite conservation laws 



In,t + Fn,x + G n , y — 0, n — 1, 2, • 



(5.18) 
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In the equation (5.18), the conversed densities I' n s are given by formula 
(5.17), the first fluxes F' n s are given by 

Fi = hihpx + (h 2 - xha)h - {xah x + 2yh i )d^ 1 h i y - Qahxd' 1 ^ + hzd^h, 
F 2 = hil 2 fix + {h 2 - xh 3 )I 2 - {xahi + 2yh i )d^ l I 2 ^ y - Ylh\hl 2 - Uahihd' 1 1 2 , y 
- Gahid^hv + 2/i 3 <9~ 1 / 2 , 

n 

F n = -6/ii^/ fc (4 + i_fc + a5 a : 1 4_fc j2/ ) - 2hi ^2 hi jh + (h 2 - xhz)I n + hJn^ 

k=l i+j+k=n 

- (xah 2 + 2yh 3 )d~ 1 I ni y + nh^d' 1 ^ - 6ahil n+ i, n = 3, 4, • • • . 
and the second fluxes G' n s are given by 

G\ = ?>hia 2 d~ l Ii^, G 2 = 3h\al 2 + ?>hiot 2 d~ l I 2: y, 

n 

G n = Shia ^2 hh-k + 2,hia 2 d~ x I n ^ n = 2, 3, • • • . 
fc=i 

The first equation of (5.18) is exactly the vc-KP equation (5.1). The 
expressions F' n s and G' n s indicate that the fluxes of the vc-KP equation 
are not local. To summarize, the vc-KP equation (5.1) is complete inte- 
grable, since it admits bilinear Backlund transformation, Lax pair and 
infinite conservation laws. 

6. General variable-coefficient KP equation 

Consider a general vc-KP equation [3D] 

(u t + hiu 3x + h 2 uu x ) x + h 3 u 2y + h A u xy + (h 5 + h 6 y)u 2x + h 7 u x = 0, (6.1) 
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where hi = hi(t), i = 1,2, ••• ,7 are arbitrary functions with respect 
to time variable t. The equation (6.1) include many special variable- 
coefficient equations in physics, such as cylindrical KdV equation [22] 

u t + uu x + u 3x + Y t u x = 0, (6.2) 



cylindrical KP equation [UJ 



1 3d 2 

(w* + + h 2 uu x ) x + — m x + -j^u 2y = 0, (6.3) 

generalized cylindrical KP equation [42], S3] 

1 3d 2 

(«t + /iiw 3x + h 2 uu x ) x + — w x + — u 2y + r(t)w xy + [f(t) + g{t)y]u 2x = 0. 

(6.4) 

Here we attempt to find the integrability condition that the equation 
(6.1) possesses bilinear representation, Backbend transformation, Lax 
pair, Darboux covariant Lax pair and infinite conservation laws. 

6.1. Bilinear representation 

By introducing a potential field q 

u = c{t)q 2x , 

with c = c(t) is free function to be determined, the resulting equation 
(6.1) for q (integrating with respect to x twice) reads 

E(q) - q xt + hq 4x + -h 2 q\ x + h 3 q 2y + h A q xy + (h 5 + h 6 y)q 2x 

* (6-5) 

+ (h 7 + d t ln c)q x = 0, 
which can be expressible as P-polynomials 

E(q) = P xt {q) + hPteiq) + h 3 P 2y (q) + h 4 P xy (q) + (h 5 + h 6 )P 2x (q) 

+ (h 7 + d t In hih 2 l )q x = 0, 

(6.6) 
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if one chooses c = 6h\h 2 1 and use the formula (2.4). By application of 
the transformation 

q = 2lnF u = cq 2x = l2hih 2 l {\wF)2 x 

and using the property (2.3), then the equation (6.6) gives the bilinear 
representation for the vc-KP equation (6.1) as follows 

Starting form this bilinear equation, we can get multi-solutions, for ex- 
ample, the regular one-soliton like solution is 

ah ^ 1/ 2 ,2 kx + s y + u 

u = bmrin k sech , 

2 2 

in which A: is a arbitrary constant, while s = s(t) and u = u(t) are two 
function with respect to £, given by 

s(t) = k J h Q dt, u(t) = - J (k 3 hi + k~ 1 s 2 h3 + sh4 + h 7 + dt\nhih 2 ~ 1 )dt. 

6.2. Backlund transformation and Lax pair 

We now consider bilinear Backlund transformation and Lax pairs of 
the equation (6.1). Let q' and q be two solutions of the equation (6.5), 
we consider the following two-field condition 

E(q') - E{q) = {q 1 - q) xt + h^q' - q) 4x + 3h 1 {q' + q) 2x (q - q) 2x + h 3 {q r - 

+ h A (q' - q) xy + (h 5 + h 6 y)(q' - q) 2x + {h 7 + d t In h x h 2 l ){q' - q) x = 0, 

(6.7) 

which may produce the required bilinear Backlund transformation under 
an appropriate additional constraint. 
On introducing two new variables 

v = (q'-q)/2, w = (q' + q)/2 (6.8) 
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we rewrite the condition (6.7) as the form 

E(q') - E(q) = v xt + h x v Ax + §h 1 v 2x w 2x + h 3 v 2y + h 4 v xy + (h b + h 6 y)v 2x 
+ (h 7 + a t In h x h 2 ^v* = d x [y t {v) + hiy 3x {v, w)} + w) = 0, 

(6.9) 

with 

R(v, w) = 3/iiWronskian[^2x(^, w),y x (v)} + h 3 v 2y + /^u^ + (h 5 + /&62/)^2. 

+ (/i 7 + d t \nhih 2 l )v x . 

In order to express R(v, w) as the x- and ^/-derivative of ^-polynomials, 
we choose the constraint 

y y (v) + a(t)y 2x (v, w) = A, (6.10) 

where a = ai(£) is to be determined. Direct calculation show that 
R(v, w) = 2>hi\v 2x - a~ l [hzW2 Xl y + (2/13 - 3/iia 2 )v x v xy + 3/iia 2 v 2 x^] 

+ ^4^xy + (/15 + hy)v2x + (/17 + In hih 2 r )v x , 

(6.11) 

which can be expressible R(v,w) as the x-derivative of a linear combi- 
nation of ^-polynomials 

R(v, w) = d x [3hi\y x - 3ahiy 3x (v, w) + h 4 y y + (h 5 + h 6 y)y x ], 

if we take a simple constraints 

hj + In h\h 2 1 = 0, /i3 = 2/2,3 — 3/iia 2 = 3/iia 2 , 

namely, 

h 2 = §h l e$ h7d \ 3a 2 = h 3 h^\ 

From (6.9)-(6.11), one infers that 
w) + - A = 0, 

WtM + w) - 3ay xy {v, w) + 3A^»] + h A y v (6.12) 

+ (h 5 + = 0, 
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which leads to bilinear Backlund transformation of variable coefficient 
KP equation 

(D 2 X + aD y - X)F • G = 0, 

[A + hi(D 3 x - 3aD x D y + 3\D X ) + /^A, + (h 5 + + (3]F ■ G = 

by using the property (2.2). 

By using the Hopf-Cole transformation v = \nip and the formulae 
(2.5) and (2.6), the system (6.12) can be linearized into a Lax pair 

(ad y + Lx)ip = (ad y + dl + q 2x - = 0, 

(d t + L 2 )i/j [d t + Ah r dl - h A a- l d 2 x + (6h iq2x + h 5 + h 6 y)d x (6.13) 

+ 3/iig 3x - Shiaq xy - h±aT x q 2x + /14a -1 A] V» = 0, 
whose integrability condition 

[ady + L h d t + L 2 ]^ = 

is satisfied if u is a solution of the vc-KP equation (6.1) and a t = 0, or 
equivalently, h^hi 1 = constant. This Lax pair can be used to construct 
Darboux transformation, inverse scattering transformation for soliton 
solutions to the vc-KP equation (6.1). 

6.3. Darboux covariant Lax pair 

Let us go back to the vc-KP equation (6.1) and the associated Lax 
pair (6.13). Assume that is a solution of the following Lax pair 

(ady + Li)(f) = A0, Li = dl + q 2x , 

(6.14) 

(d t + L 2iC0Y )(J) = 0, L 2;CO v = hidl + bid 2 x + b 2 d x + 63, 
where 61, b 2 and 63 are functions to be determined. It is shown that the 
gauge transformation 

T = <j)d x (t)^ 1 = d x — a, a = d x ln(f) (6.15) 
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map the operator ad y + L\(q) onto a similar operator: 

T{ad v + L x {q))T- 1 = ad y + L 1 (q = q + Aq) with Ag = 21n0. 

It is suffices to verify that such transformation (6.15) map the d t + L 2iC0Y 
into similar one 

T(d t + L 2jCOV )T- 1 = dt + L 2 , C ov, ^2,cov = hid 3 x + b x dl + ~b 2 d x + 6 3 , (6.16) 
where 6j, j = 1,2,3 and L 2jCO v satisfy the covariant conditions 
6, = bj(q) + A6j = bj{q + Ag), j = 1, 2, 3. 

L 2l cov{q) = L 2i cov{q = q + Ag) , Ag = 2 In 0. 
It follows from (6.15) and (3,16) that 

A&i = 61 - 61 = 0, A6 2 = 12hia x + 61 a; + A6i(J, 

(6-17) 

A6 3 = 12/iicr 2a; + 2a x b 1 + A6 2 cr + b 2 , x . 
We require 61,62 and 63 in the differential polynomial form of potential 
filed q 

bj = Fj(q, q x , q 2x , q 3x , "•), j = 1, 2 

such that 

AF 7 = Fj(q + Ag, ^ + Aq x , q 2x + Ag 2x , • • • ) - Fj(q, q x , q 2xi ■ ■ ■ ) = Abj, 

with Ag rx = 2(\nq) rx , r = 1, 2, • • • . 

From eigenvalue equation in (6.14), we get the relation 

qzx = ~aa xy - (a x + a 2 ) x , 

on account of which, solving the system (6.17) yields 
61 =ci(y,t), b 2 = 6hiq 2x + c 2 (y,t), 

(6.18) 

63 = 3hi(q 3x - aq xy ) + c x q 2x + c s (y, t), 
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with ci(y,t),C2(y,t) and c^(y,t) being arbitrary functions with respect 
to y and t. From (6.14) and (6.18), we then find a Darboux covariant 
evolution equation 

(d t + L 2 ,cov)0 = 0, 

^2,cov = + cid 2 x + {Qhiq 2x + c 2 )d x + 3h x q 3x - 3hiaq 2x + c 3 . 

In particular, if setting 

ci = -/i 4 a _1 , c 2 = /15 + /i62/ 5 c 3 = h A \, 
the Darboux covariant operator L 2 ^ C0Y reduce to operator L2, namely 

^2, cov = ^2- 

Therefore the Lax pair (6.13) is Darboux covariant under constraint 
hsh^ 1 = constant. 

In a similar way, we can get higher operators 

L p , cov (q) = h x dl + hd^ 2 + • • • + ftp, = 3, 4, • • • 

which are Darboux covariant with respect to L\ step by step, so as to 
produce higher order members of the vc-KP hierarchy. 

6.4. Infinite conservation laws 

Finally, we turn to construct the conservation laws of vc-KP equation. 
For this purpose, we expect re-decompose the two- filed condition (6.7) 
into x- and ^/-derivative of 3^-polynomials. We return to revisit R(v, w) 
in two-field condition (6.9) and write it as another form 

R(v, w) = [?>hi\v x - 3ahiv x v y + (h 5 + h 6 y)v x ] x + (h 3 v y + h A v x + 2>hiav 2 x ) 

(6.19) 
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It follows from the relations (6.9), (6.10) and (6.19) that 
y2x{v, w) + ay y (v) - A = 0, 

d t y x (v) + d x [hiy 3x (v, w) - 3hiay y (v, w) + 3hi\y x (v) + (h 5 + h 6 y)y x (v)} 
+ d y [h 3 y y (v) + 3h ia y x {v) 2 + h A y x {v)] = 0, 

(6.20) 

which is slightly different from (6.10) and can produce desired conserva- 
tion laws. Especially we don't need the constraint 3a 2 = h^h^ 1 . 
By introducing a new potential function 

and it follows from the relation (6.8) that 

v x = r/, w x = q x + r\. (6.21) 

Substituting (6.21) into (6.20) yields a coupled system 

j) x + rf + ad~\ + q 2x - s 2 = 0, (6.22) 
r\ t + d x [hi{j] 2x + Qr]£ 2 - 2jf - Qar]d~\) + (h 5 + /ie2/)»7] 

(6.23) 

+ d y (h 3 d x 1 T] y + Shiar] 2 + /1477) = 0, 
where we have used the equation (6.22) to get the equation (6.23) and 

set A = e 2 . 

Substituting the expansion 

00 

r) = £ + Y, I ^^---) £ ~ n ( 6 - 24 ) 

n=l 

into the equation (6.22), equating the coefficients for power of e, then 
we obtain the recursion relations for I n as follows 

h = -\<hx = ~\^ hid \ h = \e! hidt (u 2x + ad^uy), 

1 " , ( 6 - 25 ) 

J n = ~ g ( J n,x + hh-k + Oid x I n ,y), Tl = 2, 3, • • ■ , 
fc=l 



37 

Again substituting (6.24) into (6.23) and comparing the coefficients for 
power of e provide us infinite conservation laws 

I n ,t + F njX + G n , y = 0, n = 1, 2, • • • . (6.26) 

In the equation (6.26), the conversed densities I' n s obtained by recursion 

formulas (6.25), and the first fluxes F' n s are expressible in I' n s 

Fi = hih^x - Qahilf - Qahid' 1 ^ + (h 5 + hey)I h 

F 2 = hih,2x ~ 6hih(2I 2 + ah + ad' 1 ^) - §hiad~ l h^ + (h 5 + h 6 y)I 2 , 

n 

F n = — 6/ii ^ I k {I n +\-k + oi-d^In-i^y) — 2/ii ^ IJjh + hiI ni 2 X 

k=l i+j+k=n 

- §ahid~ l I n +i, y + {h + h 6 y)I n , n = 3, 4, • • • . 
and the second fluxes G^s are given by 

Gi = hzd^h^ + 6ahil 2 + Vi, 

G n = Shia ^2 hl n -k + 3/iia 2 ^ 1 /^, n = 2, 3, • • • . 

fc=i 

The first equation of the conservation law equation (6.26) is exactly the 
vc-KP equation (6.1). Taking the boundary condition of p into account, 
the equation (6.26) implies that I' n s, n = 1,2,- •• constitute infinite 
conserved densities of the vc-KP equation (6.1). To this end, we remark 
that as application of these results, all equations (6.2)-(6.4) are complete 
integrable under the constraint h 2 = 6hief h7dt , since they possess bilinear 
Backlund transformation, Lax pair and infinite conservation laws. 
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